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One may identify the general properties of the neutrino mass matrix by generating many 
random mass matrices and testing them against the results of the neutrino experiments. 

There are three hght, active neutrinos whose fields Ve, Vfi, and Vr are left handed; 
there also probably are three right-handed neutrinos whose fields fre, J^r/x, and v^t 
do not participate in the electroweak interactions and are said to be sterile. Fields 
that mix must transform in the same way. So it makes sense to combine the three 
left-handed (two-component) neutrinos z/e, Vfi^ and Vt with the left-handed fields 
Se, s^, and s-r that are the charge conjugates of the right-handed neutrino fields i^^e, 
Vr)i, I'rr, i-C, Sg = — i(T2 v*^ = {^1^2^ "^^rei); ^tc. The six left-haudcd neutrino fields 
form a six- vector, N = [v^, f^, 1/^, Sg, s^, s,-). The only mass term available for two 
left-handed neutrino fields Ni and Nj is M.ij {NiiNj2 — Ni2Nji) + h.c. Because the 
fields Ni anticommute, the six-by-six complex mass matrix 



is symmetric; here T means transpose. In most models the quantities Mij are the 
mean values of Higgs fields in the vacuum. Because the complex mass matrix M 
is symmetric, it may be factored Ai — U M by a matrix U that is unitary, and 
a matrix M that is diagonal and positive.^ The elements of Ad are the masses 
of the neutrinos. The vector N„i of six mass eigenfields is Nm — N . The mass 
eigenfields are the normal modes of the action. They are Majorana neutrinos. 

The six-by-six complex mass matrix M. involves 42 real parameters. One may 
explore this large space and test various properties of the mass matrix JV[ by gen- 
erating many random mass matrices that share those properties, and by computing 
the extent to which they fit the results of the neutrino experiments. For instance, 
one may define an ang le by sin^ = Tv{E^ E + F'<F)/Ti{2DW + E^E + F^F) 
and test whether Xi, can be very small. ^ This angle characterizes where the six neu- 
trinos lie on a continuum that extends from three purely Dirac neutrinos, Xi, = 0, 
to six purely Majorana neutrinos, Xi, = 7r/2. When the angle Xi, is small, as 
in a theory in which B ~ L is nearly conserved, the six masses coalesce into 
three pairs of nearly degenerate masses, and the six neutrinos form three nearly 
Dirac neutrinos, a condition that has been called pseudo-Dirac.^ This B — L or 
small-Xi/ property explains the large mixing angles and tiny mass differences seen 
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experimentally and allows the masses of the neutrinos to lie in the eV range with 
J2i 'Tii 8 eV which is a cosmological bound."* To test this property, I used the 
software package LAPACK^ to factorize 10,000 random mass matrices A4 every 
parameter of each of which was a complex number z = x + iy with x and y chosen 
randomly and uniformly from the interval [— leV, leV]. The solar neutrinos were 
taken to have an energy of 1 MeV, and the probability P{i'i — > i^i') — ~^ ^i')]"^ 

with A{vi Uii) = Y^. U*,jUij exp{—im'jL/{2E)) was averaged over one revolution 
of the Earth about the Sun. The atmospheric neutrinos were averaged over the 
atmosphere and over energies in the range of 1-30 GeV weighted by the flux of 
atmospheric muon neutrinos as a function of energy and local zenith angle given 
by the Bartol group. ^ The resulting scatter plots fit the gross features of the solar 
and atmospheric experiments quite well with sinxi, = 0.003 when two additional 
properties are added. ^ The first is a constraint on inter-generational mixing which 
I imposed by suppressing the singly off-diagonal matrix elements of D, E, and F 
by 0.2 and the doubly off-diagonal matrix elements by 0.04. The second is a weak 
mass hierarchy which I implemented by scaling the z,j-th elements of the sub- 
matrices E,F, and D of the mass matrix Ai, Eq. (|l|), by the factors f{i) * f{j) 
where / — (0.2,1,2). These three properties also lead to general agreement with 
the results of the LSND and KARMEN2 experiments. If the physical mass matrix 
Ai shares these properties, and if MiniBooNE can achieve a sensitivity of 0.001 for 
j/^ — > Vg and a precision of 0.01 for j/^ P^, then it has a good chance of seeing 
the appearance of De and the disappearance of Vfj_.^ 

Another test of the properties of the neutrino mass matrix M is provided by 
the passage of high-energy neutrinos through the Earth.'' In momentum space the 
Dirac equation for the 24- vector {N,ia2N*) of neutrino fields is 

(E + g-p + V -Mt \( N \ 
"^^y -M E-a-p-V )\ia2N* ) 

in which the hermitian matrix V represents the interaction of the neutrinos with 
the electrons and quarks of the Earth. The matrix V is diagonal with elements 
V = {Gf/V2) (2iVe - Nn, -iV„, -Nn, 0, 0, 0) in which and iV„ are the number 
densities of electrons and neutrons. By neglecting neutrinos of positive helicity, 
antineutrinos of negative helicity, and the terms V'^ and [y, A^], one finds for the 
square of the mass operator the six-by-six hermitian matrix M(p)j_ = Ad^Ai ± 
2\p\V = Md{p)^ U in which U is another unitary matrix and the plus sign is 
for neutrinos and the minus sign for anti-neutrinos. I used the PREM model^ to 
divide the Earth into 82 shells each with its approximate values of the densities 
and Nn- My code computes the amplitude for a neutrino to propagate through the 
Earth, which is the path-ordered exponential 



A = V<cxp 



I dxM{p)\/{2\p\) 



(3) 

as a product of up to 82 factors of the form C/^ exp [i Ax Mrf(p)^/(2|p|)] U. For 16 
random mass matrices M. with small x^,, little inter-generational mixing, and a weak 
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mass hierarchy, the probabihties PsKi^^p, ^ i^fj.) that a muon neutrino of energy 100 
GeV would traverse the Earth at a given azimuthal angle 9 and enter the Super 
Kamiokande detector are plotted in Fig. Because high-energy interactions with 
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Figure 1: The probability that a 100-GeV muon neutrino would traverse the Earth are 
plotted against the cosine of the azimuthal angle 6 for 16 random mass matrices M with 
sina::,^ — 0.003, little inter-generational mixing, and a weak mass hierarchy. 

matter suppress active-sterile neutrino oscillations, the probability Psk(^'p ^ i^p.) 
remains close to unity in 13 of the cases. But in three cases it drops quite far below 
unity. Further work is required to determine what additional properties, if any, 
would cause such random mass matrices to display the oscillations and neutral- 
current events observed by the Super-Kamiokande collaboration.^ 
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